function answer = polynom_i1(x, y, z)
 answer=0.011370335673695078216443743576154*x - 0.1107647591873588307950661516333*y + 0.055584370666661293839209044780425*z - 0.029583623327723371711716195543019*x*z + 0.0069442777792777698854664044831247*x*(x*z - (57*z)/10 + 1/5) - 0.015782449498358567921514555643465*z*((61*y)/25 + (61*z)/25) - 0.0059167246655446743423432391086038;
 end
